Driven linear modes: Analytical solutions for finite discrete systems 
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We have obtained exact analytical expressions in closed form, for the linear modes excited in finite 
and discrete systems that are driven by a spatially homogeneous alternating field. Those modes are 
extended for frequencies within the linear frequency band while they are either end-localized or 
end-avoided for frequencies outside the linear frequency band. The analytical solutions are resonant 
at particular frequencies, which compose the frequency dispersion relation of the finite system. 
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Introduction.- The calculation of the linear modes 
which can be excited in finite and discrete systems that 
are driven by a spatially homogeneous alternating (AC) 
field is of great interest. That problem arises in the low- 
amplitude limit of several model equations of physical 
systems, which are comprised of elements that are either 
directly or indirectly coupled (i.e., the so called magne- 
toinductive systems). In that limit, it is a common prac- 
tice to omit the nonlinear terms ending up with a linear 
non-autonomous system of equations that accept linear 
wave solutions called in the following driven linear modes 
(DLM). Such a linearized problem may arise, for exam- 
ple, from the AC driven Frenkel-Kontorova model^— , 
which has been widely used to model rf-driven paral- 
lel arrays of Josephson junctions (JJs)^~— , or from more 
general models of AC driven anharmonic lattices with 
realistic potentials^. It may also arise from models of 
inductively coupled AC driven rf-SQUID arrays^—, from 
models of inductively coupled intrinsic Josephson junc- 
tions supplied by AC current^—, and also from nonlin- 
ear magnetic metamaterial models that are comprised of 
split-ring resonators placed in an AC magnetic fielcU^— . 
In the latter case, those linear waves are known as (lin- 
ear) magnetoinductive wave s 19 ' 20 , and they have phonon- 
like dispersion curves 2 ^ and many prospects for device 
applications 2 ^—. In the present work we present analyt- 
ical solutions that we have obtained for DLMs in finite 
and discrete systems, with frequencies either in the linear 
wave band (LWB) obtained in the usuall way, or outside 
that band. We find that the former ones are extended 
modes as it is expected, while the latter ones are either 
end-localized or end-avoided modes. The results are il- 
lustrated using as a paradigmatic example the driven FK 
chain 2 , which in a standard normalization form reads 
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where qi is the ith 'coordinate' whose interpretations de- 
pends on the particular system to which Eq. (fTJ) is re- 
lated, a is the loss coefficient, C is the coupling con- 
stant, fac is the amplitude of the sinusoidal driver with 



frequency w, and i = 1,...,N, with N being the total 
number of elements in the chain. 

The linearized problem.- Linearization of Eq. (fTJ) gives 

q\ + qi = C(qi-x - 2q { + q i+1 ) + f ac shr(wi), (2) 

where the losses have been omitted for simplicity. The 
earlier equation accepts solutions of the form qi{t) — 
Qi sin(wt), where Qi is the time-independent mode am- 
plitude at site i. By substitution of qi{t) into Eq. ([2]) we 
obtain a system of stationary equations 
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where s and k are real parameters given by 
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In general, s and k depend on the parameters of the 
particular model. Since we consider finite systems, Eq. 
Q should be implemented with open-ended boundary 
conditions, i.e., Qo = Qn+i — 0, to account for the 
termination of the structure at both ends. The formal 
solution of Eq. ([3]) can be written as 



Q = kST 1 !, 



(5) 



where Q and 1 are TV— dimensional vectors with com- 
ponets Qi and 1, respectively, and S _1 is the inverse 
of the N x N coupling matrix S. The latter is a real, 
symmetric tridiagonal matrix that has diagonal elements 
equal to unity, while all the other non-zero elements are 
equal to s. The need to find the inverse of tridiagonal 
matrices like S arises in many scientific and engineering 
applications. Recently, Huang and McColl 2 ^ related the 
inversion of a general trigiagonal matrix to second or- 
der linear recurrences, and they provided a set of very 
simple analytical formulae for the elements of the inverse 
matrix. Those formulae lead immediately to closed forms 
for certain, relatively simple trigiagonal matrices, like S. 

Analytical solutions.- Following Huang and McColl 26 
we can show that the elements of S _1 are given by [i > j) 
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for -1/2 < s < 1/2, where 



for < s < +1/2, 
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and 



sin(jfl') sin[(A - i + 1)0'] 

M ~ n, ZIZTTTt , ^/./i > ( 8 ) 



'ij ^ sin0' sin[(JV + l)0'] ' 
for s > 1/2 or s < —1/2, where 



The elements of S _1 for i < j are given by the same ex- 
pressions as in Eqs. ([5J and (ITU1) by simply interchanging 
the indices i and j (due to the symmetry). Although the 
above expressions are valid for any N, in the following 
we assume that N is even. 

For obtaining the solution Q, we first write Eq. (|3]) as 
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It turns out that the sum in the right-hand-side of the 
earlier equation can be be calculated in closed form, using 
in an appropriate way the summation formulae 2 ^ 
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for -1/2 < s < +1/2, and the formulae 2 -^ 
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for s > 1/2 or s < —1/2. Then, using several identi- 
ties for hyperbolic and trigonometric functions, we arive, 
after long and tedious calculations to the solution 
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for -1/2 < s < 
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for s > +1/2, and 



where 



= sin — cos 



(24) 



for s < -1/2. 

The case where s = ±1/2 should be treated separately. 
However, the solution of Eqs. (|3|) for that specific value 
of the coupling parameter can be obtained by calculating 
the limit s — > 1/2 of the expressions Eqs. (| 1 5 [) - (|TH]) and 
Eqs. P0 > ([2T ]) (for s = +1/2). We find that the limits of 
those expressions are the the same, so that Qi(s = 1/2) 
is given by 



Qi—odd ^ 
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In the same way, by calculating the limit s — > —1/2 of 
the expressions Eq. (JT5J) and Eq. ((33J) (for s = -1/2) 
gives for the Qj(s = —1/2) the expression 



Qi = ni(N -i + 1). 



(27) 



At this point we have given the solutions of Eqs. ([3]) for 
all real values of the coupling parameter s. 

Illustrative examples and discussion.- For illustration 
we show in Figs. 1 and 2 (for s > and s < 0, respec- 
tively) plots of Qi as a function of i for an FK chain with 
N = 50, C = 0.5, f ac = 0.02, and several frequencies. 
We should note that those results where checked numer- 
ically. The linear dispersion relation (LDR) for Eq. @ 
;*2 
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FIG. 1: (color online) Driven linear modes for the finite 
Frenkel-Kontorova chain for C = 0.5, f ac = 0.02, N = 50, 
and s > 0. The continuous lines serve as a guide to the eye. 
(a) End-localized modes for u> = 2.646 (s = 0.1, k = -0.004; 
black circles), uj = 1.915 (s = 0.3, ft = —0.012; red squares), 
to = 1.744 (s = 0.48, ft = -0.0192; green triangles), (b) Spe- 
cial case of ui = u> max (s = +0.5, ft = —0.02). (c) Extended 
modes for u> = 1.717 (s = 0.52589, k = -0.021; black circles), 
lo = 1.459 (s = 3.8781, k = -0.0155; red squares). 
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FIG. 2: (color online) Driven linear modes for the finite 
Frenkel-Kontorova chain for C = 0.5, f ac = 0.02, N = 50, 
and s < 0. The continuous lines serve as a guide to the 
eye. (a) End-avoided modes for ui — 0.2773 (s = —0.26, 
k = 0.0104; black circles), uj = 0.5774 (s = -0.3, ft = 0.012; 
red squares), ui = 0.9428 (s = —0.45, ft = 0.018; green trian- 
gles), (b) Special case of uj = uj m i n (s — —0.5, ft = 0.02). (c) 
Extended modes for u = 1.016 (s = -0.51668, ft = -0.02067; 
black circles), u = 1.089 (s = -0.61493, ft = -0.02457; red 
squares) . 



where k is the normalized wavenumber. The extends 
from w min = 1 to io max = y/l + AC, having gaps be- 
low and above uj m in and io max , repsectively. It is easy 
to see that the frequencies within the LWB correspond 
to either s < —1/2 or s > +1/2, for u) m i n < uj < uj s 
and uj s < uj < uj max , respectively, where ui s = yl + 2C. 
Several DLMs are shown in Figs, la, lb, and lc, cor- 
responding to < s < +1/2, s = +1/2 and s > +1/2, 
respectively. We observe that the modes in Fig. la, with 
frequencies above the LWB, exhibit weak localization at 
the end-points of the chain. The mode shown in Fig. lb 
correspond to frequency uj = Lu max , exactly on the upper 
bound of the LWB, while those shown in Fig. lc are ex- 
tended stationary solutions, with frequencies in the LWB. 
The DLMs shown in Figs. 2a, 2b, and 2c, correspond to 
> s > — 1/2, s = -1/2 and s < -1/2, respectively. 
The modes shown in Fig. 2a, with frequencies below the 
LWB, exhibit end-avoidance. The mode shown in Fig. 
2b correspond to frequency ui = Lu m i m , exactly on the 
lower bound of the LWB, while those shown in Fig. 2c 
are extended stationary solutions. From Figs. 1 and 2 we 
observe that the amplitude of the DLMs with frequencies 
in the LWB increases with decreasing frequency. 



Frequency dispersion and concluding remarks.- Note 
that those amplitudes are determined uniquelly from the 
parameters of the system. Even though we have chosen 
a rather weak driving amplitude f ac , the DLMs close to 
the lower bound of the LWB attain rather large ampli- 
tude (> 1 in some cases in normalized units). Moreover, 
from the analytical solutions we infer that the DLMs are 
resonant for particular values of frequencies in the LWB. 
The frequencies of those resonances can be obtained by 
zeroing the denominators of the solutions for the DLMs 
in the Eqs. (17]), ([19]), (22]), and (24]). Specifically, by 
zeroing the denominators ZX£ and DTp , for s > +1/2 and 
s < —1/2, respectively, we get for the resonant values of 
s the relations 



= < 2 cos 



2mn 



(29) 



for s > +1/2, where the integer < m < (n + l)/4, while 



-2 cos 



(2m + 1)tt 
n + 1 
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for s < —1/2, where the integer < m < (n — l)/4. The 
total number of resonances, which provide the frequency 
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FIG. 3: (color online) Frequency dispersion (eigen-frequencies 
vs. mode number m) for the finite Frenkel-Kontorova chain 
with N = 20, and C — 0.5. The black curve gives the fre- 
quency dispersion for the infinite system while the red cir- 
cles correspond to the eigen-frequencies of the periodic sys- 
tem with N = 20. The green squares and the blue rhombii 
are the eigen-frequencies of the finite system with open-ended 
boundary conditions for m even and odd, respectively. 

dispersion for a finite chain, are n/2. The corresponding 
resonance frequencies uj m are obtained from the first of 
Eqs. Q, as 

u m = yJu 2 s + (C/s*), (31) 

with given by either by Eq. ([H} (for s > +1/2) or by 
Eq. (|3U)) (for s < —1/2). Since we ignored the loss term 
in Eq. the mode amplitudes go to infinity at the 

resonant frequencies. We should note that both in Figs. 



lc and 2c, the values of to were chosen in between two 
neighboring resonances. If losses were taken into account, 
the mode amplitudes would still reach considerably high 
values, even for very weak driver. That example ques- 
tions the validity of any linear approximation for DLMs 
at some part of the LWB, which for the particular model 
lies between oj s and uj m i n . That situation get worse and 
worse as the lower bound is approached, and the large 
amplitude modes are expected to become modulationally 
unstable for moderately high driving amplitudes. 

We have presented analytical solutions in closed form 
for the DLMs in finite and discrete systems whose ele- 
ments are coupled with their nearest neighbors. It should 
be stressed that these solutions are generally valid for any 
such system with the appropriate forms s and k on its pa- 
rameters. The presented examples, using the driven FK 
model show that the DLMs are extended for frequencies 
in the LWB and end-localized (end-avoided) for frequen- 
cies above (below) the LWB. Those particular features 
of DLMs are due to the choise of the open-ended bound- 
ary conditions, which modifies the local potential that 
an element feels close to the end-points. The resonance 
frequencies w m , obtained from Eq. (|3"TT) for all relevant 
values of m, provide the frequency dispersion for the fi- 
nite chain, which differs slightly from that of the periodic 
system (Fig. 3). The difference increases with decreasing 
number of elements. Preliminary calculations show that 
the localized states shown here are stable and, moreover, 
they may evolve into dissipative surface states when the 
nonlinearity sets in and the losses are included^. 
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